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Sunyaev-Zel’dovich Effect Using the Skewness of the CMB Temperature Distribution 
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We present a detection of the unnormalized skewness ( T"(h )) induced by the thermal Sunyaev- 
Zerdovich (tSZ) effect in filtered Atacama Cosmology Telescope (ACT) 148 GHz cosmic microwave 
background temperature maps. Contamination due to infrared and radio sources is minimized by 
template subtraction of resolved sources and by constructing a mask using outlying values in the 
218 GHz (tSZ-nuil) ACT maps. We measure (T 3 (n)) = —31 ± 6 pK 3 (measurement error onlj ) 
or ±14 pK 3 (including cosmic variance error) in the filtered ACT data, a 5a detection. We show 
that the skewness is a sensitive probe of as, and use analytic calculations and tSZ simulations to 
obtain cosmological constraints from this measurement. From this signal alone we infer a value of 
as = 0.78io!o4 (68% C.L.) to'.il (95% C.L.). Our results demonstrate that measurements of non- 
Gaussianity can be a useful method for characterizing the tSZ effect and extracting the underlying 
cosmological information. 


I. INTRODUCTION 

Current observations of the cosmic microwave back- 
ground (CMB) anisotropies on areminute scales using 
experiments such as the Atacama Cosmology Telescope 
[ACT; 1-3] and the South Pole Telescope [SPT; 4, 5] 
probe not only the primordial microwave background 
fluctuations sourced 13.7 billion years ago, but also 
measure secondary anisotropies caused by more recent 
and less distant physical processes. Such secondary 
anisotropies are induced by infrared (IR) dustv galax- 
ies and radio sources, gravitational lensing, and the 
Sunyaev-Zerdovich (SZ) effect. The SZ effect [6, 7] arises 
due to the inverse Compton scattering of CMB photons 
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off high energy electrons located predominantly in hot 
gas in galaxy clusters (the intrar-cluster medium, or ICM). 
This scattering modifies the spectrum of CMB photons in 
the direction of a cluster in a way that depends on both 
the thermal energy contained in the ICM (the thermal 
SZ effect) as well as the peculiar velocity of the cluster 
with respect to the CMB rest frame (the kinetic SZ ef- 
fect). The kinetic SZ effect simply increases or decreases 
the amplitude of the CMB spectrum in the direction of 
a cluster, but the thermal SZ (tSZ) effect modifies the 
CMB spectrum in a frequency-dependent manner. The 
tSZ effect is characterized by a decrease (increase) in the 
observed CMB temperature at frequencies below (above) 
218 GHz in the direction of a galaxy cluster due to in- 
verse Compton scattering. The thermal effect is generally 
at least an order of magnitude larger than the kinetic 
effect for a typical massive cluster at 148 GHz. Mea- 
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surements of the tSZ signal, which is proportional to the 
integrated ICM pressure along the line of sight, can be 
used zo observe the high redshift universe, constrain cos- 
mological parameters - in particular erg, the variance of 
matter fluctuations on scales of 8 Mpc jh - and probe 
baryonic physics in the ICM. 

The tSZ signal has so far primarily been studied either 
by directly resolving individual clusters in arcminute- 
scale CMB maps [8-12] or by measuring it statistically 
through its presence in the small-scale CMB power spec- 
trum [13, 14]. However, it is by no means obvious that 
the power spectrum is the best way to characterize the 
statistical properties of the tSZ field. Indeed, measuring 
the tSZ signal in the power spectrum is challenging be- 
cause there are many other sources of CMB power on ar- 
cminute scales: primordial CMB fluctuations, CMB lens- 
ing, instrumental noise, dusty star- forming IR galaxies, 
and radio sources. In order to disentangle these contri- 
butions to the power spectrum and isolate the amplitude 
of the tSZ signal, a sophisticated multifrequency analysis 
is required, which involves modeling the power spectrum 
contribution of each of these components in at least two 
frequency bands. 

In this paper we instead measure the tSZ signal using 
the unnormalized skewness of the filtered temperature 
fluctuation (T 3 (n)). This quantity has the significant 
advantage that, unlike measurements of the tSZ effect 
through the power spectrum, its measurement does not 
require the subtraction of Gaussian contributions, be- 
cause it is only sensitive to non-Gaussian signals with 
non-zero skewness. The primordial CMB (which is as- 
sumed to be Gaussian on these scales) and instrumen- 
tal noise (which is Gaussian) hence do not contribute 
to it. In addition, CMB lensing and the kinetic SZ ef- 
fect do not induce skewness (as they are equally likely 
to produce positive and negative fluctuations), and so 
do not contribute either. The primary contributions 
to this quantity are thus only the tSZ effect and point 
sources. These signals have a different frequency depen- 
dence. Furthermore, the tSZ signal contributes nega- 
tive skewness, whereas radio and IR point sources con- 
tribute positive skewness. These characteristics allow the 
tSZ signal to be effectively isolated and studied, as first 
pointed out in [15]. 

Measurements of the skewness also possess significant 
advantages from an astrophysical perspective. A con- 
sistent problem plaguing studies of the tSZ power spec- 
trum has been theoretical uncertainty in the ICM elec- 
tron pressure profile [16-18], especially in the low-mass, 
high-redshift groups and clusters that contribute much 
of the signal. As discussed in the following section in de- 
tail, the tSZ skewness signal is dominated by character- 
istically higher-mass, lower-red shift clusters than those 
that source the power spectrum signal. The ICM astro- 
physics for these objects is better constrained by X-ray 
observations and they are less sensitive to energy input 
from non-gravitation al sources [18, 19]. Thus, the theo- 
retical systematic uncertainty in modeling the tSZ skew- 


ness is correspondingly lower as well. In addition, at 148 
GHz, dusty star-forming galaxies are less present in 
massive, low-redshift clusters (which contribute more to 
the skewness) than in high-redshift groups and clusters 
(which contribute more to the tSZ power spectrum) [20]. 
Thus, we expect the correlation between tSZ signal and 
dusty galaxy emission, which can complicate analyses of 
the tSZ effect, to be smaller for a measurement of the 
skewness. 

Moreover, the tSZ skewness scales with a higher power 
of ag than the tSZ power spectrum amplitude. This re- 
sult is precisely what one would expect if the signal were 
dominated by higher-mass, rarer objects, as the high- 
mass tail of the mass function is particularly sensitive 
to a change in erg. This provides the prospect of tight 
constraints on cosmological parameters from the skew- 
ness that are competitive with constraints from the power 
spectrum. 

In this paper, we first explain the usefulness of the 
skewness as a cosmological probe by theoretically deriv- 
ing its scaling with as as well as the characteristic masses 
of the objects sourcing the signal. Subsequent sections of 
the paper describe how we measured this skewness in the 
ACT data. We describe how the ACT temperature maps 
are processed in order to make a reliable measurement of 
the unnormalized skewness due to the tSZ effect, and dis- 
cuss how contamination from IR dusty galaxies and radio 
point sources is minimized. We report the measurement 
results and discuss how the errors are calculated. Finally, 
we discuss the cosmological constraints and associated 
uncertainties derived from this measurement. 

We assume a flat ACDM cosmology throughout, with 
parameters set to their WMAP5 values [21] unless other- 
wise specified. All masses are quoted in units of M q //i, 
where h = i?o/(100kms -1 Mpc -1 ) and Hq is the Hubble 
parameter today. 

II. SKEWNESS OF THE TSZ EFFECT 

In this section, we investigate the N th moments of the 
pixel probability densitv function, (T n ) = (T(n) N ), fo- 
cusing on the specific case of the unnormalized skewness 
(T 3 ). We show that the unnormalized skewness (T 3 ) has 
a steeper scaling with as than the power spectrum ampli- 
tude and is dominated by characteristically higher-mass, 
lower-redshift clusters, for which the ICM astrophysics 
is better constrained and modeled. As explained earlier, 
these characteristics make tSZ skewness measurements a 
useful cosmological probe. 

In order to calculate the N th moment of the tSZ field, 
we assume the distribution of clusters on the sky can be 
adequately described by a Poisson distribution (and that 
contributions due to clustering and overlapping sources 
are negligible [-22]). The N th moment is then given by 

a) 
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where T(0,M,z) is the tSZ temperature decrement at 
position 6 on the sky with respect to the center of a 
cluster of mass M at redshift z: 

T'(0;M,z) — g(v)TcMB~—^s 
m e c z 

x J P^yp + ^izm^M.z^dl, (2) 

where g(v) is the spectral function of the tSZ effect, 
cIa(z) is the angular diameter distance to redshift z, and 
the integral is taken over the electron pressure profile 
P e (r; M, z) along the line of sight. 

For a given cosmology, Eqs. (1) and (2) show that there 
are two ingredients needed to calculate the N th tSZ mo- 
ment (in addition to the comoving volume per steradian 
dV/dz , which can be calculated easily): (1) the halo mass 
function dn(M, z)/dM and (2) the electron pressure pro- 
file P e (r; M, z) for halos of mass M at redshift z. We use 
the halo mass function of [23] with the redshift-dependent 
parameters given in their Eqs. (5)-(8). While uncertain- 
ties in tSZ calculations due to the mass function are often 
neglected, they may be more important for the skew- 
ness than the power spectrum, as the skewness is more 
sensitive to the high-mass exponential tail of the mass 
function. We estimate the uncertainty arising from the 
mass function by performing alternate calculations with 
the mass function of [24], which predicts more massive 
clusters at low redshift than [23] for the same cosmol- 
ogy. As an example, the predicted skewness calculated 
using the pressure profile of [16] with the mass function 
of [23] is » 35% lower than the equivalent result using 
the mass function of [24]. However, the derived scalings 
of the variance and skewness with as computed using [24] 
are identical to those found below using [23]. Thus, the 
scalings calculated below are robust to uncertainties in 
the mass function, and we use them later to interpret 
our skewness measurement. However, we rely on cosmo- 
logical simulations to obtain predicted values of the tSZ 
skewness. We do not consider alternate mass functions 
any further in our analytic calculations. 

We consider three different pressure profiles from [16, 
17, 25] in order to evaluate the theoretical uncertainty in 
the scaling of the tSZ skewness with as- These profiles 
differ in how they are derived and in the ICM physics 
they assume. They thus provide a measure of the scatter 
in the scalings of the variance and skewness with as due 
to uncertainties in the gas physics. 

Finally, in order to make a faithful comparison be- 
tween the theory and data, we convolve Eq. (2) with 
the Fourier-space filter described in the subsequent data 
analysis sections of this paper. In addition, we imple- 
ment a high-mass cutoff in Eq. (l) ’that roughly corre- 
sponds to the pixel fluctuation cutoff used in the data 
analysis. These steps cannot be neglected, as the filter 
and cutoff reduce the predicted tSZ skewness amplitude 
bv up to 95% compared to the pure theoretical value [26] . 
Most of this reduction is due to the filter, which modestly 
suppresses the temperature decrement profile of a typi- 


cal cluster; this suppression strongly affects the skewness 
because it is a cubic statistic. 

The analytic theory described above determines the 
scaling of the N th tSZ moment with cr 8 . In particular, 
we compute Eq. (1) with N = 2 and N = 3 for each 
of the chosen pressure profiles while varying as- The 
scalings of both the > r ariance ( N = 2) and the skewness 
(N — 3) with as are well-described by power laws for 
each of these profiles, (T 2 ’ 3 ) oc a % 2,3 . For the profile of 

[16] , we find <*2 = 7.7 and 0:3 = 10.8; for the profile of 

[17] , we find a 2 = 7.9 and <23 = 11.0; and for the profile of 
[25], we find 0:2 — 7.5 and 0:3 = 10.5. Note that the scal- 
ing of the variance matches the scaling of the tSZ power 
spectrum amplitude that has been found by a number 
of other studies, as expected (e.g., [25, 27]). The scaling 
of the unnormalized skewness is similar to that found by 
[28], who obtained a 3 = 10.25. Also, note that the skew- 
ness scaling is modified slightly from its pure theoretical 
value [26] due to the Fourier-space filter and high-mass 
cutoff mentioned above. The overall conclusion is that 
the skewness scales with a higher power of ag than the 
variance (or power spectrum). We use this scaling to 
derive a constraint on ag from our measurement of the 
skewness below. 

In addition, we compare the characteristic mass scale 
responsible for the tSZ skewness and tSZ variance (or 
power spectrum) signals. Analytic calculations show 
that the tSZ power spectrum amplitude typically receives 
« 50% of its value from halos with M < 2-3 x 10 14 Mq//i, 
while the tSZ skewness receives only « 20% of its am- 
plitude from these less massive objects. This indicates 
that the clusters responsible for the tSZ skewness signal 
are better theoretically modeled than those responsible 
for much of the tSZ power spectrum, both because mas- 
sive clusters have been observed more thoroughly, and 
because more massive clusters are dominated by gravita- 
tional heating and are less sensitive to non-linear energy 
input from active galactic nuclei, turbulence, and other 
mechanisms [18, 19]. We verify this claim when inter- 
preting the skewness measurement below, finding that 
the systematic theoretical uncertainty (as derived from 
simulations) is slightly smaller than the statistical error 
from the measurement, though still non- negligible. 

III. MAP PROCESSING 
A. Filtering the Maps 

The Atacama Cosmology Telescope [1-3] is a 6m tele- 
scope in the Atacama Desert of Chile, which operated 
at 148, 218, and 277 GHz using three 1024-element ar- 
rays of superconducting bolometers. The maps used in 
this analysis were made over three years of observation 
in the equatorial region during 2008-2010 at 148 GHz, 
and consist of six 3° x 18° patches of sky at a noise level 
of « 21 /iK arcmin. In our source mask construction we 
also use maps of the same area made in 2008 at a fre- 
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quency of 218 GHz. The maps were calibrated as in [29]. 
We apodize the maps by multiplying them with a mask 
that smoothly increases from zero to unity over 0.1° from 
the edge of the maps. 

Although atmospheric noise is removed in the map- 
making process, we implement an additional filter in 
Fourier space to remove signal at multipoles below t = 
500 {t is the magnitude of the Fourier variable conjugate 
to sky angle). In addition, we remove a stripe for which 
-100 < ^dec < 100 along the Fourier axis correspond- 
ing to declination to avoid contamination bv scan noise. 
Furthermore, to increase the tSZ signal-to-noise, we ap- 
ply a Wiener filter which downweights scales at which the 
tSZ signal is subdominant. This (non-optimal) filter is 
constructed by dividing the best-fit tSZ power spectrum 
from [13] by the total average power spectrum measured 
in the data maps, i.e. C\ sz /C\ ot . For multipoles above 
£ = 6 x 10 3 , the tSZ signal is completely dominated by 
detector noise and point sources, and hence we remove 
all power above this multipole in the temperature maps. 
The final Fourier-space filter, shown in Fig. 1, is normal- 
ized such that its maximum value is unity. As it is con- 
structed using the binned power spectrum of the data, it 
is not perfectly smooth; however, we apply the same fil- 
ter consistently to data, simulations, and analytic theory, 
and thus any details of the filter do not bias the interpre- 
tation of our result. After filtering, the edges of the maps 
are cut off to reduce any edge effects that might occur 
upon Fourier transforming despite apodization. Simula- 
tions verify that no additional skewness is introduced by 
edge effects into a trimmed map. 

B. Removing Point Sources 

In order to obtain a skewness signal due only to the tSZ 
effect, any contamination of the signal by point sources 
must be minimized. These objects consist of IR dusty 
galaxies and radio sources. We use two approaches to 
eliminate the point source contribution: template sub- 
traction and masking using the 218 GHz channel. 

In the template subtraction method [30], which we 
use to remove resolved point sources (mainly bright ra- 
dio sources), sources with a signal-to-noise (S/N) greater 
than five are first identified in a match-filtered map. A 
template with the shape of the ACT beam is then scaled 
to the appropriate peak brightness of each source, and 
this profile is subtracted from the raw data. The pro- 
cess is iterated following the CLEAN algorithm [31] until 
no more sources can be identified. We verify that this 
procedure does not introduce skewness into the maps 
(e.g., through oversubtraction) by checking that similar 
results are obtained using a different procedure for reduc- 
ing source contamination, in which we mask and in-paint 
pixels which contain bright sources with S/N > 5 [32]. 

We take a second step to suppress the lower-flux, unre- 
solved point sources (mainly dusty galaxies) that remain 
undetected by the template subtraction algorithm. At 
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FIG. 1. The Wiener filter applied to the ACT temperature 
maps before calculating the unnormalized skewness. This fil- 
ter upweights scales on which the tSZ signal is large compared 
to other sources of anisotropy. 


218 GHz, dusty galaxies are significantly brighter than 
at 148 GHz and the tSZ effect is negligible. We construct 
a dusty galaxy mask by setting all pixels (which are ap- 
proximately 0.25 arcmin 2 ) to zero that have a tempera- 
ture in the 218 GHz maps larger than a specified cutoff 
value. This cutoff is chosen to be 3.2 times the standard 
deviation of the pixel values in the filtered 148 GHz map 
(3.2cr). This procedure ensures regions with high flux 
from dusty galaxies are masked. We also set to zero all 
pixels for which the temperature is lower than the nega- 
tive of this cutoff, so that the masking procedure does not 
introduce spurious skewness into the lensed CMB distri- 
bution, which is assumed to have zero intrinsic skewness. 
The mask is then applied to the 148 GHz map to re- 
duce the point source contribution. Simulations ([33] for 
IR sources) verify that the masking procedure does not 
introduce spurious skewness into the 148 GHz maps. 

Finally, all pixels more than twelve standard deviations 
(12 a) from the mean are also removed from the 148 GHz 
maps. Due to the ringing around very positive or negar 
tive pixels caused by the Wiener filter, the surrounding 
eight arcminutes of these points are also masked. This 
additional step slightly increases the S/N of the skew- 
ness measurement by reducing the dependence on large 
outliers, enhances the information content of low mo- 
ments by truncating the tail of the pixel probability den- 
sity function, and ensures that any anomalous outlying 
points from possiblv mis-subtracted bright radio sources 
do not contribute to the skewness signal. Overall, 14.5% 
of the 148 GHz map is removed by the masking proce- 
dure, though the removed points are random with respect 
to the tSZ field and should not change the signal. 
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FIG. 2. Histogram of the pixel temperature values in the FIG. 3. Likelihood of the skewness measurement described in 
filtered, masked ACT CMB temperature maps. A Gaussian the text (with measurement errors), 
curve is overlaid in red. 


IV. RESULTS 
A. Evaluating the Skewness 

We compute the unnormalized skewness of the filtered 
and processed 148 GHz maps by simply cubing and av- 
eraging the pixel values in real space. The result is 
(X 3 ) = —31 i 6 /iK 3 , a 5 a deviation from the null result 
expected for a signal without any non-Gaussian compo- 
nents. The skewness of the CMB temperature distribu- 
tion in our filtered, processed maps is visible in the pixel 
value histogram shown in Fig. 2 (along with a Gaussian 
curve St for comparison). The likelihood corresponding 
to our measurement of the skewness is shown in Fig. 3. 

The measurement error on the skewness includes only 
Gaussian sources of noise. We calculate this error by 
using map simulations that consist of Gaussian random 
fields with the same power spectrum as that observed in 
the data, including noise and beam effects. This esti- 
mate thus does not include the error resulting from cos- 
mic variance (i.e. sample variance) of the non-Gaussian 
tSZ signal. 

We calculate errors that include cosmic variance by 
constructing more realistic simulations. To construct 
such simulations, we add maps with simulated tSZ signal 
from [34], which assume = 0.8, to realizations of a 
Gaussian random field which has a spectrum such that 
the power spectrum of the combined map matches that 
observed in the ACT temperature data. Given the sim- 
ulated sky area, we obtain 39 simulated maps, each of 
size 148 deg 2 . By measuring the scatter amongst these 
patches and appropriately scaling the error to match the 
237 deg 2 of unmasked sky in the processed map, we ob- 
tain a full error (including cosmic variance) on the un- 
nornu lized skewness of 14 /iK 3 . The scatter of skew- 
ness values measured from each of the simulated maps 
is consistent with a Gaussian distribution. The estimate 
for the full error is coincidentally the same as the stan- 


dard error, 14 /iK 3 , estimated from the six patches into 
which the data are divided. The full error is used below 
in deriving cosmological constraints from the skewness 
measurement. 


B. The Origin of the Signal 

Is the skewness dominated by massive clusters with 
large tSZ decrements - as suggested by theoretical con- 
siderations described earlier - or by more numerous, less 
massive clusters? To investigate this question, we mask 
clusters in our data which were found in the 148 GHz 
maps using a matched filter as in [35]. All clusters de- 
tected above a threshold significance value are masked; 
we vary this threshold and measure the remaining skew- 
ness in order to determine the origin of the signal. 

Fig. 4 shows a plot of the signal against the cluster 
detection significance cutoff. We include calculations us- 
ing both the full cluster candidate catalog obtained via 
matched filtering as well a catalog containing only clus- 
ters confirmed optically using the methodology of [36] 
on the SDSS Stripe 82 [37]. The SDSS Stripe 82 imag- 
ing data cover « 80% of the total map area, and thus 
some skewness signal will necessarily arise from objects 
not accounted for in this catalog. The results for these 
two catalogs agree when masking clusters with S/N > 7, 
but differ slightly when masking lower S/N clusters. This 
effect is likely due to the small shortfall in optical follow- 
up area as well as a small number of false detections (i.e., 
impurity) in the candidate clusters that have not yet been 
optically followed up. 

Using either catalog, Fig. 4 implies that just under 
half of the tSZ skewness is obtained from clusters that 
lie below a 5 a cluster detection significance, while the 
remainder comes from the brightest and most massive 
clusters. The results of [35] suggest that clusters de- 
tected at 5 a significance are roughly characterized by a 
mass M 500 = 5 x 10 14 M© /A, where M 500 is the mass 
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Minimum S/N of masked detected clusters 

FIG. 4. Plot of the skewness signal as a function of the mini- 
mum S/N of the clusters that are masked (this indicates how 
many known clusters are left in the data, unmasked). The 
blue line is calculated using the full cluster candidate catalog 
obtained via matched filtering, while the green line uses a cat- 
alog containing only optically-confirmed clusters [37], Both 
lines have identical errors, but we only plot them for the green 
line for clarity. Confirmed clusters source approximately two- 
thirds of the signal, which provides strong evidence that it is 
due to the tSZ effect. Note that one expects a positive bias 
of ~ 4 fxK for the S/N = 4 point of the blue line due to im- 
purities in the full candidate catalog masking the tail of the 
Gaussian distribution. 


enclosed within a radius such that the mean enclosed 
density is 500 times the critical density at the cluster red- 
shift. This value corresponds to a virial mass of roughly 
M = 9 x 10 14 M©//i, which was also found to be the 
mass detection threshold for high-significance ACT clus- 
ters in [36], Fig. 4 thus demonstrates that roughly half 
of the tSZ skewness signal is due to massive clusters with 
M > 10 15 M 0 /fr. The theoretical calculations described 
earlier give similar results for the fraction of the signal 
coming from clusters above and below this mass scale, 
which is significantly higher than the characteristic mass 
scale responsible for the tSZ power spectrum signal. 

Finally, the positive value in the full, candidate cat- 
alog line shown in Fig. 4 when masking clusters above 
S/N = 4 is consistent with zero. When masking at this 
level (with the candidate catalog which contains some 
impurities), we slightly cut into the negative pixel val- 
ues in the Gaussian component of Fig. 2, leading to a 
small spurious positive skewness. For the points we plot, 
we calculate that this bias is only non-negligible for the 
S/N = 4 cut, where it is « 4 fiK. This bias effectively 
explains the small positive offset seen in Fig. 4. How- 
ever, we discuss positive skewness due to any possible 
residual point source contamination below. Overall, the 
dependence of the measured skewness on cluster masking 
shown in Fig. 4 provides strong evidence that it is caused 
by the tSZ effect. 



FIG. 5. A test for Ift source contamination: similar to the 
blue line in Fig. 4, but with a range of values of the cutoff 
used to construct an IR source mask in the 218 GHz band. 
Any cutoff below « 3.2a gives similarly negative results and 
thus appears sufficient for point source removal, where a — 
10.3 fiK is the standard deviation of the 148 GHz maps. For 
comparison, the standard deviation of the 218 GHz maps is 
» 2.2 times larger. The percentages of the map which are 
removed for the masking levels shown, from the least to the 
most strict cut, are 0.7, 2.5, 8.4, 14.5, 23.7 and 36.6%. 


C. Testing for Systematic Infrared Source 
Contamination 

Despite our efforts to remove point sources, a small 
residual point source contamination of the signal could 
remain, leading to an underestimate of the amplitude 
of the tSZ skewness. To investigate this systematic er- 
ror source, we vary the level at which point sources are 
masked in the 218 GHz maps (the original level is 3.2 
times the standard deviation of the pixel values in the 
filtered 148 GHz map (3.2<r), as described above). The 
results of this test are shown in Fig. 5, which uses the 
full catalog of cluster candidates as described in Fig. 4, 
since the optical catalog does not yet cover the entire 
ACT map. Note that at « 3.2a the signal asymptotes 
to the skewness value quoted earlier and no longer be- 
comes more negative with stricter masking levels. This 
result suggests that masking at the 3.2a level sufficiently 
removes any contamination bv IR sources, and stricter 
cuts will just reduce map area unnecessarily. 

To further estimate the residual point source contami- 
nation in the 148 GHz maps, we process simulations of IR 
sources from [33] (with source amplitudes scaled down by 
1.7 to match recent observations) with the same masking 
procedure as that applied to the data, creating a mask in 
a simulated 218 GHz map, and applying it to a simulated 
IR source signal at 148 GHz. We find a residual signal 
of (T 3 ) = 4.7 /iK 3 . However, this is an upper bound 
on the contamination, as the template subtraction pro- 
cedure should substantially reduce the spurious skewness 
due to IR sources. 

We also investigate a linear combination of the 148 and 
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218 GHz maps that should have minimal IR source levels, 
namely, an appropriately scaled 218 GHz map subtracted 
from a 148 GHz map. Assuming that the spatial distribu- 
tion of the point sources is not affected by the difference 
in observation frequency between 148 and 218 GHz and 
a single spectral index can be applied to all sources, a 
simple factor of w 3.2 [13] relates a point source’s signal 
in the two different frequency bands. We find that the 
appropriate linear combination (subtracting 1/3.2 times 
the 220 GHz map from the 148 map) produces a signal 
in agreement with that resulting from the previously de- 
scribed masking procedure, although the additional noise 
present in the 218 GHz maps slightly reduces the signif- 
icance of the detection. 


V, COSMOLOGICAL INTERPRETATION 

To obtain cosmological information from the measured 
amplitude of the unnormalized skewness, we compare 
our results with two different sets of tSZ simulations 
[33, 34]. Both sets of simulations are run with <r 8 = 0.8, 
but differ in their treatment of the ICM. The simula- 
tion cf [34] is a fully hydrodynamic cosmological simu- 
lation that includes sub-grid prescriptions for feedback 
from active galactic nuclei, star formation, and radiative 
cooling. The simulation also captures non- thermal pres- 
sure support due to turbulence and other effects, which 
significantly alters the ICM pressure profile. The simu- 
lation of [33] is a large dark matter-only TV-body simu- 
lation that is post-processed to include gas according to 
a polytropic prescription. This simulation also accounts 
for non-thermal pressure support (though with a smaller 
amount than [34]), and matches the low-redshiffe X-ray 
data presented in [17]. 

We perform the same filtering and masking as that ap- 
plied to the data in order to analyze the simulation maps. 
For both simulations, the filtering reduces the signal by 
« 95% compared to the unfiltered value. For the simula- 
tions of [34], we measure ( T 3 ) = -37 /aK 3 , with negligi- 
ble errors. However, this value is complicated by the fact 
that these simulations only include halos below 2 = 1 . 
An analytic estimate for the skewness contribution due 
to halos with z > l from Eq. (1) gives a 6% correction, 
which yields (T 3 ) = —39 /aK 3 . For the simulations of 
[33], we measure (T 3 ) = — 50 /aK 3 , with errors also neg- 
ligible for the purposes of cosmological constraints. 

In order to compare our tSZ skewness measurement 
with these simulation results, we use the scaling of the 
skewness with <r 8 calculated analytically above, denoted 
by a 3 .. The analytic calculations give a 3 = 10.8,11.0, 
and 10.5 for the pressure profiles of [16], [17], and [25], 
respectively, while [28] found a 3 = 10.25. For simplicity, 
we assume a scaling of a 3 = 10.5. The resulting con- 
straint on as is fairly insensitive to the exact value of 
this scaling for 10 < a 3 <12. 

We estimate the value of cr 8 from our skewness mea- 


surement as simply 


-as = of 




-| 1/OE3 


Ufa) 5 ] 


( 3 ) 


where quantities with the superscript S and D cor- 
respond to simulation and data, respectively. Using 
the simulation values in combination with the skew- 
ness measurement and full error (including cosmic vari- 
ance) described earlier, we compute Eq. (3) and obtain 
a B = 0.78±g;jg (68% C.L.) (95% C.L.) for [34] and 

ag = 0.76+°;°3 (68% C.L.) (95% C.L.) for [33]. For 

comparison, the final release from the Chandra Cluster 
Cosmology Project found a& = 0.803 dr 0.0105, assum- 
ing = 0.25 (there is a strong degeneracy between 
cr 8 and for X-ray cluster measurements that probe 
the mass function) [38] . Perhaps more directly com- 
parable, recent studies of the tSZ power spectrum have 
found as = 0.77 ± 0.04 (statistical error only) [13] and 
as = 0.807d=0.016 (statistical error and approximately es- 
timated systematic error due to theoretical uncertainty) 
[39]. 

To evaluate the theoretical systematic uncertainty in 
the amplitude of the filtered skewness due to unknown 
ICM astrophysics, we test the effect of different gas 
prescriptions by analyzing simulations from [34] with all 
forms of feedback, radiative cooling and star formation 
switched off, leading to an adiabatic ICM gas model. 
For these adiabatic simulations we find ( T 3 ) = —56 fiK 3 
(after applying the 6% correction mentioned earlier), 
which for the skewness we measure in our data would 
imply cr 8 = 0.76t°;° 3 (68% CL)jft?jj (95% CL). Turning 
off feedback and all sub-grid physics is a rather extreme 
case, so the systematic theoretical uncertainty for a 
typical simulation with some form of feedback should 
be slightly smaller than the statistical error from the 
measurement, though still non-negligible. This .contrasts 
with measurements of cr 8 via the tSZ power spectrum, 
for which the theoretical systematic uncertainty is 
comparable to or greater than the statistical uncertainty 
[13, 39]. As highlighted earlier, this difference can 
be traced to the dependence of the power spectrum 
amplitude on the ICM astrophysics within low-mass, 
high-redshift clusters. The skewness, on the other hand, 
is dominated by more massive, lower-redshift clusters 
that are less affected by uncertain non-gravitational 
feedback mechanisms and are more constrained by 
observations. Nonetheless, as the statistical uncertainty 
decreases on future measurements of the tSZ skewness, 
the theoretical systematic error will quickly become 
comparable, and thus additional study of the ICM 
electron pressure profile will be very useful. 
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VI. CONCLUSIONS 


As the thermal Sunyaev-ZePdovich field is highly non- 
Gaussian, measurements of non-Gaussian signatures such 
as the skewness can provide cosmological constraints that 
are competitive with power spectrum measurements. We 
have presented a first measurement of the unnormalized 
skewness (T 3 (n)) in ACT CMB maps filtered for high 
signal to noise. As this is a purely non-Gaussian signa- 
ture, primordial CMB and instrumental noise cannot be 
confused with or bias the signal, unlike measurements 
of the tSZ power spectrum. We measure the skewness 
at 5 <t significance: (T 3 (n)) = —31 ± 6 /xK 3 (measure- 
ment error) . Including cosmic variance increases the er- 
ror to zb 14 fi K 3 . Using analytic calculations and sim- 
ulations to translate this measurement into constraints 
on cosmological parameters, we find <r 8 = 0.78io!o4 (68% 
C.L.) (95% C.L.), with a slightly smaller but non- 

negligible systematic error due to theoretical uncertainty 
in the ICM astrophysics. This detection represents the 
first realization of a new, independent method to mea- 
sure (7s based on the tSZ skewness, which has different 
systematic errors than several other common methods. 
With larger maps and lower noise, tSZ skewness mea- 
surements promise significantly tighter cosmological con- 


straints in the near future. 
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